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ABSTRACT 


In the present work transfer finite element equations are 
developed for a four noded rectangular element and an eight 
noded isoparametric e 1 emen t . The s e formulations are next used for 
solving two dimensional and ax i-symme t r i c thermal and stress 
analysis problems employing transfer FEM. A rectangular domain and 
a pressure vessel with hemi-s pher i cal heads is considered for two 
dimensional and ax i-symme t r i c problems respectively. Results are 
obtained for different geometrical and natural boundary 
in both two dimensional and ax i -symme t r i c problenis. 
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CHAPTER - 1 


INTRODUCTION 

Exact methods fail to provide a solution for most of the 
real life engineering problems. This is more true when 
components involved are non-uniform. Approximate methods have to 
be used for obtaining reasonably accurate solutions for such 
problems. Such methods suffer from the limitations of 
satisfying the geometrical and natural boundary conditions right 
in the beginning. To overcome these limitations numerical methods 
like finite difference, transfer matrix method and finite element 
methods have been developed. In the present work the latter two 
are relevant and are discussed here. 

Transfer matrix method was the first one to be developed. 
Its earliest application was to the torsional vibrations of 
shafts by Holzer Cn in WZI. In 1944, after twenty three years 
it was extended to beam vibrations by Myklestad C23 . It was 

further developed by many research workers including Pestel and 

Leckie L31 who wrote the text C33 in 1963. Large references are 

given in this text C33 ; and problems like vibrations of 

Timoshenko beams, circular arches, turbine generator shafts, 
shear panels etc. are studied. Method is also applied to static 
problems of beams, rotating disks and stability problems of 
beams. It may be noted that ninety nine percent of this text is 
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devoted to one dimensional problems. They also tried to extend 
this transfer matrix approach to plate vibrations by replacing 
the plate by an equivalent network of beams based on the 
Hernnikoff model. This type of modelling is very artificial and 
questionable. A more general procedure for getting transfer 
matrices for plate and shell problems is given by Pestel C4I]. 

Next text, Pilkey and Chang C5I1 , dealing with this transfer 
matrix came after fifteen years i.e. in 1978. Probably the last 
and latest to come is Tesar and Fillo C6I] in 1988. These texts 
deal with plate problems more extensively. Thus it may be noted 
that probably only four texts have come in this twentieth century 
so far. 

These transfer matrices are usually obtained from statics or 
by using results that are tabulated in engineering hand books. 
They are also obtained from the closed form solutions of the 
governing differential equations of the problem, but quite a 
difficult and many a times impossible job. In addition many 
analytical methods like Laplace transformation, exponential 
expansion, matrizant, product integral C3,5D and numerical 
methods C5,6Z1 are used to obtain the transfer matrices. 

Since fifties when modern electronic computers became 
available, finite element method got developed extensively. 
Now-a-days it is the most widely accepted technique for the real 
life enginering problems. Probably more than three hundred texts 
have been written since 1967. 

Beauty of these methods lie in the fact that the boundary 
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conditions are introduced in the end. In transfer matrix methods 
one multiplies transfer matrices one after the another, so this 
method does not tax the memory of the computers. But as said 
earlier, accurate transfer matrices are not that easy to obtain. 
In finite element method, finite element matrices are added and 
the size of assembled matrices goes on increasing with the 
increase in number of elements. Thus one faces serious 
limitations of computer memory when using this method for more 
involved problems. 

Thus in view of the above mentioned limitations of transfer 
matrix method and finite element methods, a new approach was 
developed by Dokainish C7J in 1972 for vibrations of rectangular 
plates using triangular elements with three nodes (incompatible 
elements). This approach combines the advantages of both these 
methods i.e. derivation of finite element matrices is very 
straight forward (these are changed to transfer matrices in a 
very simple way) and very small computer memory is needed. This 
new approach is called transfer finite element method (TFEM) in 
this work. It may be pointed out that correct FEM equations were 
not used in this work but correct TFEM equations were obtained 
using the knowledge of statics. McDaniel and Eversole CS3 
proposed a similar approach in 1977 for a stiffened plate 
structure. This method was probably improved by Pilkey C9I in 
Technical Report in 1979. Sankar and Hoa C10T in 1980 offered an 
approach in which an extended transfer matrix relating the state 
vectors which consists of state variables displacements and 
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forces and their derivatives with respect to frequency were used. 

This approach was also independently developed by Gupta C11I1 
in late 82. Correct FEM equations were used to develop transfer 
matrices. This method was developed for free vibrations of 
one-dimensional problems of axial rods, Euler beams, Timoshenko 
beams (Single fourth order differential equation and two 
simultaneous second order coupled differential equations), 
rotating Euler and Timoshenko beams and two dimensional problems 
of rectangular isotropic and anisotropic membranes using 
rectangular four noded elements. This was also developed for 
static problems of one dimensional axial rods and Euler beams 
and two dimensional steady state heat conduction problem in a 
rectangular region using triangular elements with three nodes. 
Results obtained by classical approach and TFEli were found to be 
identical for all practical purposes. 

Tesar C1EI] applied this method to non-linear vibrations of 
shell structures in 1983, te chn i cal report in Slovak, and extended 
this work C13I] in 1985. Ogha and Shigematsu C14I] used it in 1983 
for bending and buckling problems and developed a technique 
for treating the structure with intermediate conditions. 

TFEM was used by Subramanya C153 for flexibility analysis of 
pipes in 1985. Pipes having branches and intermediate 
constraints were taken care of. Various loadings were 

considered. This method was used by Goel C163 for finding 
critical speeds of rotors in 1986 and also by Mathew C17D for 
studying flexural vibrations of rotating shafts in 1992. Mathew 
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C173 studied flexural vibrations of damped shafts including the 
effect of shear deformation, rotary inertia and gyroscopic 
actions using four second order differential equations. Critical 
speeds, complex natural frequencies and unbalance responses were 
obtained. All types of geometrical and natural boundary 
conditions, including bearing cross coupled properties were 
incorporated exactly. 

Gupta C18G used this method in 1987 for vibration analysis 
of cable networks and beam grillages. Orthogonal and 
non~orthogonal cable networks were considered. Transverse 
vibrations of beam grillages in bending and combined-torsional 
modes with various' types of boundary conditions were analysed. 

This method has been developed considerably by Tesar [180 and 
applied to many linear and nonlinear complex problems like thin 
walled beams, corrugated and sandwich panels, shells, 
elastic-plastic behaviour, aeroelastic response and so on. Tesar 
C8I] calls this method as finite element transfer matrix (FETM). 
It may be noted that here also correct FEM equations are not used 
in this text C63 but correct TFEM equations are obtaied using the 
knowledge of statics later. In two dimensional problems only 
four no.ded elements are considered. 

In TFEM, the domain is divided into number of strips. Strip 
directions are so chosen that minimum number of elements are 
needed in these strips. Each strip is divided into number of 
finite elements. Governing finite element matrices are assembled 
for this strip in the usual way. The boudary conditions of top 
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and bottom edges are incorporated. These equations are 
rearranged to obtain the transfer matrix for this strip. These 
transfer matrices are multiplied, boundary conditions of the left 
and the right ends are applied, and thus obtaining the first 
state vector. Other state vectors are obtained in the usual way. 

In the present work TFEM equations have been developed for 
two-dimensional heat conduction and stress analysis problems 
using four noded rectangular and eight noded isoparametric 
elements. Convective heat transfer boundary conditions have also 
been considered. All this is given in Chapter 2. 

Above is also done for axi-symmetr i c heat conduction and 
stress analysis problems and is given in Chapter 3. Here a case 
of unequal adjacent strips has also been considered. Method is 
also extended to transient heat conduction. 

To validate the above, results are obtained for different 
problems and are given in Chapter 4. 

Chapter 5 concludes the present work. 
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CHAPTER - 2 

STEADY STATE THERMAL AND STRESS ANALYSIS OF TWO DIMENSIONAL 

PROBLEMS 


In this chapter TFEM eii^uations have been developed frorri 
FEM equations for two dimensional steady state heat conduction 
(with various boundary conditions) and stress analysis problems 
using four noded rectangular and eight noded isoparametric 
e 1 emen t s . C on ve c t i V e boundary conditions have also been 
considered. Transfer matrices have been developed for a typical 
strip(first strip)of the plane lamina. In the last transfer finite 
element method (TFEM) has beem illustrated for both thermal and 
stress analysis problems. 


2.1 Governincf dlff’erent/lal equablon^ heat, conduction 

The governing differential equation for heat conduction in a 
plane with no heat source or sink is C203 , Fig. 2.1.1 


S 

dx 


(k 


dT, 


X di 





0 


( 2 . 1 . 1 ) 



Fig. 2.1.1 Plane lamina with domain D and boundary B 

The boundary conditions can be of the types 
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o r 


T 


T 

P 


( :< , y ) 


1 

k - 5 — n 

i< Sx ;< 


I 

k - 3 — n 
y Sy y 


q + h(T - T^) = 


( 2 . 1 . 2 ) 


2.2 FEM equations^ heat conduction 
2.2.1 Without convection 

Here the field variable, temperature, is approximated over 


the element as, 
,(e) 


T 


= a+bK+cy+dxy+ 


L N^J ^ 


Tg J- = lNJ -CT} 


( ne ) 


( 2 . 2 . 1 ) 


Substituting this in equation (2.1.1), residue is obtained 


as. 


„(e) ^ £_ .. . ±_ 

^ dx ^x dx dy ^y dy 


( 2 . 2 . 2 ) 


Minimizing this residue by the Galerkin method, one gets 


JJ ^ [k 


(k . (k., dy = 0 


■ X d X 


Sy y Sy 


A 


( e ) • 


Doing intergration by parts, one gets 


(2.2.3) 


dT 


( e ) 




N. k^ - 5 - ) dx dy 

i.y y ^y 

^^(e) ar^®^ 

(2.2.4) 
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These equations in the matrix form are 


JJ (k^fN,x> LN,x J + k <N,y> LN,yJ) d:< dy 

= ^ fN.> (-q) dB 

and in short form is 


(2.2.5) 


r-, -, ( e ) r_,(ne) ^ ^(ne) 

Ck^3 -CTj- = -CqJ 


( 2 . 2 . 6 ) 


where, 


(e) 


is thermal stiffness matrix of the element 


fT>(ne) is elemental temperature vector 

is flux vector along the edges of the element. 

2. 2. 2 With convection 

Denoting temperature difference, T-T^ , by B, the governing 


differential equation (2.1.1) can be written as 

|_ (k 1^) . (k = 0 

dx X By y By 

and boundary conditions as 
<9 = & 


(2.2.7) 


o r 


k +k n + q + hd=0 

X ^x X y By y 


( 2 . 2 . 8 ) 


( e ) 

& over the finite element is taken as 
= ’ = LNJ 


(2.2.9) 


Substituting it in equation (2.2.7) one gets 


(e) 


£_ (i- ii'®’) + <k 

Bx ^^'x ^x «y y By 


( 2 . 2 . 10 ) 


Minimizing this residue by Galerkin method. 
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n [li 


<k 


. i- <k “"’)]d>< dy = 0 


'A Sx 


Sy y Sy 


< 2 . 2 . 11 ) 


Doing integration by parts, one gets 




(e) 


ite 


(e) 


rr <N. k ^ + N. k ^ ) dx dy 

JJ i,x X ^x i,y y Sy 


< e ) 


= <frN.(k 5^ n +k ^ n)dB 
J 1 K X y y 


^ N. (-h (9^®^ - q) dB 


= ^ - h N^LNJ dB Ce} ^ q N^dB (2.2.12) 


Boundary condition <2.2.8) has been invoked in equation <2.2.12). 
These equations in the matrix form are 


<k^^ -CNyX } LN,xJ + k^ -CNryJ LN,yJ)dx dy+ ^ h CN} LNJdBj -Cd} 

= ^ - q {N> dB 


< ne ) 
<2.2. 13) 


and in short form is 


t h 


o r 




{q> 


< ne ) 


<2.2.14) 


<e) 


It may be noted that matrix Ck^^l is meaningful only along 


the edges having convection 
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2. 3 TFEM equa'bions for heat conduction 

@ © (D ® 3 

2 

I I I 1 '^ 1 

Fig. 2. 3. lCa3 Rectangular lamina Fig. 2. 3. I<lb2> Strips divided 

divided Into strips Into elements 

First the domain is divided into number of finite strips as 
shown in Fig. 2.3.1(a). Strip directions are so chosen that 
minimum number of elements are needed in these strips. Strip 
could have any number of elements, but here each strip is shown 
to have two elements only, as this reduces the typing work 
considerably and enables the method to be explained very easily. 
Here each finite element is taken as four noded rectangular 
e 1 ement . 

2.3.1 FEM equa-tlons ot a strip Cwlthout convection! 

On assembling the elemental stiffness matrices and the 
flux vectors i.e. equation (2.2.6), one gets FEM equations for 
the first strip as 


© 6 © (D © 


A 




A 

5 




4 
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(2.3. 1 ) 


Transfer FEM equations are to be developed from these 
equations. So the integrals qlj^^, changed to 

*^^^65' '^^^36' this RHS of equation (2.3.1) becomes 



14 

— 


“'ll I 21 

- 


“'^^132 

+ 

"‘136 


+ 

"‘I 54 

^ I 54 

4. 

"‘l65 

■^^65 

+ 

"‘136 


(2.3.2) 


2.3.2 FEM «quAt,lons of a sst.rlp Cwl-Lh oonvactloiO 

Proceeding as in 2.3.1, the FEM equations for the first 
strip consisting of elements 1 and 2 (using equations (2.2.14)) 
are 


12 



^11 ^12 ^16 


■ ■ 



^21 ^22 ^26 


^2 


-qllg^ - 

^31 ^32 ^36 


^3 


-ql|32 + ql|3^ 


< 



> 



^4 


-ql|.i4 + ^11154 

« « m 




■^1|54 + qiUs 

^61 ^62 ^66 -I 


- ^6 - 




(E.3.3) 


2.3.3 FEM equations of a strip after applying boundary 
conditions on its top and bottom 

For developing TFEM equations of the strip, boundary 
conditions on its top and bottom edges should be applied first. 
Following four cases with different boundary conditions on top 
and bottom edges are illustrated. 

(a) Top and bottom edges insulated. 

(b) Top edge with prescribed temperature and bottom edge 
insulated. 

(c) Top and bottom edges with prescribed temperature. 

(d) Top edge exposed to fluid (convection) and bottom edge 
with prescribed temperature. 
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Case CaD Top and bottom edges insulated 


^ Q e Q (D ® 


jrr/jLi A 

A 

r £ / 1 


■/■/■ZZZ 
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5 

V7j/jr 

rr-TJ'Fj 

hrirr 


4 


Fig. 2. 3. 2: Rectangular lamina with top 


and boltom 


Insulated 




As edges 14 and 36 are insulated, equations (E.3.1) become 


"61 "62 



^<12 



^22 

^26 

m 






It 
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- 

- 
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‘•••II ■ 


'66 



■ "^1 ■ 

■ 


^2 


1 

■^3 

(i 


^4 






- ^6 - 



<1*121 

-1*132 - 1*l21 
-I'IsE 

1*l54 

1*l65 ■" 1*l54 

11 1 65 




(2.3.4) 
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Cas^Cb^ pi'^scri.b^&^l on t/op odgo sind ins^oJl j&'tod on "Lbo 

bottom ^dge 

, © - ® Tp=T* 


1 


A 




A 

5 

4 




Fig. 2. 3. 3 Rectangular lamina with bottom edge Insulated and 
top edge having prescribed ten^perature 

Applying boundary conditions of edge 36 (i.e. and 

) and edge 14 (i.e. qlj^^ =0), FEM equations for this 


strip become 


11 

*^12 

‘‘^14 

*‘•15 ^'^ 13 ^ 3 '^*‘^ 16 ''’ 6 ^ 

21 

*"'22 

*‘■24 

*‘^25 <'‘^ 23 ^^*^ 26 ■^^ 

31 

'"•32 

*^34 

45 ^ W ^*^ 36 ''^■ 

41 

*‘‘42 

*^^44 

'< 45 ^^^ 43 ^^*‘^ 46 ■'^ 

51 

^2 

K54 

‘< 55 < S 3 ■'^ S 6 ■^^ 

61 

*^’62 

*^64 

'< 65 ^‘‘^ 63 ''^'‘^ 66 -' I > . 



6:<5 



■ ■'l ■ 



■^2 


< 

■^4 

■ ^5 

> = <( 


1 


-■ 

5x1 

L 


-ql| 


21 


qilsi ~ 

q 1 j 32 

•^^32 

■^^36 


54 


■^^54 

ql |65 ■" 

q 1 1 26 

6 x 1 



( 2 . 3 . 5 ) 


It will be seen later that TFEM equations for this strip 
can be obtained from first, second, fourth and fifth equation of 
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(2.3.5) and third or sixth equation can later be used to find 
flux along edge 36. Thus writing the above mentioned four 
equations of (2.3.5) 


‘''11 

*'^12 
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14 

‘<15 

'‘<13■'^"16^^■ 


■ ■ 


^^21 

*<22 

‘‘■24 
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'‘‘^23^^‘<26■‘■^ 


^2 


‘<41 

•<42 

*<44 
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•<51 

‘<52 

■<54 

•<55 

<‘‘33''^‘<56‘^2> 


T's 








1 



‘^M21 " "^^132 
■^^154 

qiUs ■^^154 


(2.3.6) 


Case Cc!) Temperature prescribed on bottom and top edges 

3 

2 

1 

Fig. 2.3.4 Rectangular lamina with tenperature prescribed on 
opposite edges 
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Applying faoundatry conditions of edge 36 (i.e. T 

2 

‘'■6 " '•' 6 ^ ^"^36 14 (i.e. = T* and = 

equations (2.3.1) and (2.3.2) for this strip become 


^12 *^15 ^‘^11^1'^’^ 13^3'^ '"•14^4'^ ^16^6^ 



1<:22 ^33 < Kh , ) 


■ ■^ 2 ' 



*'■32 *^35 ^''•31^l'^‘‘-33^3'^^'34^4'^‘^36^6^ 




-ql|32+ql|3^ 

*'•42 *^45 ^^41^l‘^’'-43^3‘^*^44'^4'^'^46^6^ 

< 


> < 

-‘>'I-,4-^'’M54 

S 2 S 5 <S1T^S3■'^S4■'^'<56■^^ 


1 


■’*154 ■"■’Mas 

*^62 *^65 ^‘^61^l'^‘''63^3'^'^64^4'*‘'^66^6^ 


■’Mas ■"■^‘136 

L, 

* 

• M 


(2.3.7) 

It will be seen later that TFEM equations for this strip can 
be obtained from second and fifth equation of ( 2 . 3 . 7 ) and 
remaining equations can later be used to find flux across the 
edges 14 and 36. Thus writing the above mentioned two equations 
of (2.3.7) 


T 3 and 
) , FEM 



<iM54+'’M65 


(2.3.8) 
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CaseCd^ Fltiid on top edge and prescribed ten^erature on the 


bottom edge 

3 

2 

1 

Fig. 2. 3. 5 Rectangular lamina with fltiid on top edge and 


0 6® 



prescribed temperature on the bottome edge 


Applying the boundary conditions of edge “14 (i.e. = T^-T^ 
and = Tji “ noting the convection on edge 36 FEM 
equations (2.3.3) for this strip become 


^12 

^■13 


^22 

^23 

^25 

^32 

^33 

^35 

^42 

^43 

^45 

^52 

^53 

Ss 

^62 

^63 

^65 


^•16 ^^ 11 ® 1 '^^ 14®4 
^26 <^21^^^24^4 
^36^^3-1®l'^^34^4^ 
^46^^41*^1‘*'^44^4^ 
^ 56 ^ ^ 5 - 1® -1 "^^ 54 ^ 4 ^ 
^ 66 ^ ^ 61^1 "^^ 64 ^ 4 ^ 



- qiU4 - qi|2i 

4II2I ~ '^^132 

"■^^32 

I .^4 + '^^154 

’^^54 "■ ‘^^Us 

^ I 65 


(2.3.9) 


It will be seen later that TFEM equations for this strip can 
be obtained from second, third, fifth and sixth equation of 
(2.3.9), and first or fourth equation can later be used to find 
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flux along edge 14. 
equations of (2.3.9) 


Thus writing 


the 


above mentioned four 


^22 

^23 

^ 2 S 

^26 

^^ 21 ^^^ 24 ^^■ 


' ^2 ■ 

r 

^32 

^33 

^35 

^36 

<^ 31 ^^^ 34 ^^ 


^3 


^52 

^53 




i 

^5 

> < 

^62 

^63 

^65 

^66 



i 








1 



4 II 21 


(2.3.10) 


2.3.4 Transfer finite element equations 

The transfer FEM equations are now developed for all cases 
discussed in 2.3.3. It involves partitioning and rearranging of 
matrices of the strip. 


2. 3. 4. 1 TFEM equations for case da]) 

Denoting the flux ~ '^^l21 ”‘^^131 

qi|54 by qi| the FEM equations (2.3.4) are rewritten and 

partitioned as indicated. 


*''11 *^12 *^13 :‘'-14 *^15 ^■•16 

u - k - It - "If k - P 

21 22 23 .24 25 26 

• 

*^31 *^32 *^33 :*^34 *^35 ''•36 

> - - - - 

■ ■'l 

^2 

^3 . 

i • • • « • .< 

■ - qMsi ■ 

qi | 3 i 
^ 1 32 

m 

t<41 ■‘'•42 ''^43 :‘'‘^44 ''■45 *'^46 

*^51 *^52 HS ''^55 *"56 

■ ‘‘•61 ‘‘^62 ‘‘^63 :‘‘^64 ‘‘^65 ‘‘•66 - 


^4 

■^5 

- ^6 - 


ql |34 

■^^^64 

- 

( 2 . 3 . 11 ) 


where T^ , T^r T^ are temperatures and "'^^lai’’ ”^^^32 
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are 


heat fluxes for the left side of this strip; and T . , , T, 

^3 o 

temperatures and qlj^^, heat fluxes for the 

right side of this strip. Equation (2.3.11) are written in short 
form as 


CA3 

CCl 


CB3 

CD3 


CT>.^ 1 r -Cql 


fT>, 


(2.3.12) 


fq>, 


Subscript 1 outside the matrix denotes left side of the 
strip and 2 denotes right side of the strip. Equations (2.3.12) 
can be written in the following form. 


CA3 ■CT> 


CBD CT>, 


^q>.. 


CC3 •CT>^ + CDD (:T>2 


= ■Cq>, 


Equation (2.3.13) gives. 


•CT>, 


= CB3 


■^^-•Cq>.^ - CA3 •CT>.^ j 


(2.3.13) 

(2.3.14) 

(2.3.15) 


Suubst i tut ing this value of in equation (2.3.14) one 

gets , 

- CDD i:BD“''cADj ijy^ - LDl i;BD“'^^q>^ = Cql^ (2.3.16) 
Combining equations (2.3.15) and (2.3.16) in matrix form. 


one gets. 


•CT>2 



> = 

•Cq>2 


. 

m 


-CBH ^ CA3 
(CC3-CDnCB3“^CAD) 


-CBH 

-CD3CB3 



(2.3.17) 


o r 


f "■'^2 

[ 


•cq>2 

y = CTFM3.^V 

>■ 

<q>^ 


(2.3.18) 
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These are TFEM equations. Knowing tempe ratures and fluxes 
on one side of the strip, one can obtain temperatures and fluxes 
on the other side. 

2. 3. 4. 2 Transfer FEM equations for cases CbD^ Cc3 and CdD 

FEM equations for all these cases are of same nature. 
So the TFEM equations here are developed for case (b) and their 
development for other cases will be obvious. 

Rewriting equation (2.3.6), and denoting “TII 32 ~ 

-ql|3.1 and ’^^i64 

partitioned as indicated. 


■ •'^11 

S 2 

‘‘•14 *^15 



S 1 


21 

'^•21 

S2 

‘‘•24 ‘‘■25 

'S3■‘‘^S6^^ 

m m mi 

s 

^ = 4 

- ■’^31 

A. 







r * 


•<41 


‘‘•44 ‘‘^45 



s 


'’M54 

Si 


‘'•^54 

^‘‘•53S'^‘‘^56S^ 


s 








1 




(2.3.19) 

T , Tg are temperatures and -qlj^.^i' heat fluxes 

for the left side of the strip; and T^, T^ are temperatures and 
ql|^^ and are heat fluxes for the right side of the strip. 

Equations (2.3.19) are partitioned as indicated. Writing these 
in short form as, 

CAT CBT CO 

CDT CET CFl 



Subscripts 1 and 2 refer to left and right sides of the 



strip as before. 


Equation 

(2.3.20) 

can 

be 

written 

as : 


CAT 


+ CBT 

<TJ2 

+ 

CO = 

-cq>.i 

(2.3.21 ) 

and CDT 


+ CET 

TT>2 

"f 

CF} = 

^q>2 

(2.3. 22 ) 


From equation (2.3.21) one gets 

- CA3 -CT} ^ - foj (2.3.23) 

Substituting this value of CTJg in equation (2.3.22), one 


gets , 


= ^CD3-CE3CB3 ‘^CAnj CTJ^ + ^-CEnCBn fq}.^ + I^CF} - 


CEDCB3 




(2.3.24) 


Combining equations (2.3.23) and (2.3.24), one gets 



(2.3.25) 

In this equations, state vectors on left and right side are 
not of the same type. To make these of same type, these 
equations are modified as follows. 
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On combining equations (2.3.25) and (2.3.26), one gets 


■ fT>2 


-CBl ^CA3 

1 

n 

w 

u 

1 

-CB3 ’^CO 


■ ^T>; 

■cq>2 


(CD3-CE3CB3“'^CA3) 

-CE3CB3 ^ 

(fF}-CE3CB3~'^ CC> ) 

< 


1 


LOJ 

LOJ 

1 


1 


(2.3.26) 
o r 


fT>2 

< 

^ = CTFMD.^ < 

•cq>2 

1 



CT> 




> 


(2.3.27) 


These are TFEM equations, one can obtain temperatures and 
fluxes on the right side of the strip, knowing these values on 
the left hand side. 

2.4 TFEM equations for conduction with eight noded elements 

As in section 2.3, each strip is supposed to have two 
elements. But here finite elements are eight noded 
isoparametric elements (Fig. 2.4.1). The approach of section 
2.3 cannot be used directly to develop TFEM equations for this 
strip because of the intermediate nodes like 11, 12 and 13. 
These nodes are condensed out in order to develop the transfer 
finite element equations. 

For explaining this method, TFEM equations are developed for 
two cases. In the first case there is convection on one edge and 
opposite edge insulated. And in the second case there is 
convection on one edge and opposite edge has prescribed 
temperature. Extension to other boundary conditions on the 
opposite edges will be obvious. 
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Case CaD : Convection on one edge and insulat^ion on ilie other 
edge 

5 
4 

w .. 3 

2 

1 

Fig. 2.4.1 : Plane lamina with convection on one edge and 

insulation on the other edge 

Writing FEM equations for the strip using equation 
(2.2.14) . 




and qlj^Q 3 will be changed to ql|g^r <1^1 ^0,8 Also 


the boundary conditions on the top edge 5-13-10 (convection) and 
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bottom edge 1-11-6 6 ~ applied, 

equations (E.4.1) become 


^11^12 ' 
^21 ^22 


'15 


'25 


Si Se ■ 

^61 ^62 


'55 
, t 


^16‘ ■ •''^1,10 
^26 '*^2,10 


65 


^ 10,1 ^ 10,2 ^ 10,5 

^ 11,1 ^ 11 , 2 ’ ■ '^ 11,5 


^13,1 ^13,2* 


13,5 


56 

t 


.t 


5,10 


66 '‘^6,10 


^ 10 , 6 * ‘^ 10,10 
^11,6*-''^11,10 


^13;6' "^13,10 


^1,11 ••■''^1,13 
^2, 11 ■ • ’^2, 13 


^5, 11 ■■ -^5, 13 
^6,11'“'’^6,13 


^10,11 • '^10, 13 
^11,1 V^1 1,13 


'13,11 • ■ ^13, 13 



e 




& 


e 


10 


11 


13 


(2.4.2) 


Thus 
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Temperatures at the intermediate nodes of the strip are not 
required in the TFEM equation of the strip. They can be 
condensed from equations (2.4.2) C19I1. 

Equations (2.4.2) are partitioned as indicated and are 
rewritten as 


CA1T 

CA23 

CA3T 


■CT>.i ■ 



CA4T 

CAST 

CA6T 

< 

fT>2 

> = < 

I 

•Cq>2 I 

CA7T 

CAST 

CA9T 




fO> j 


Third equation of equations (2.4.3) gives. 


= CA93 


-1 


-CA73CT>.^ - CA83 fT 




(2.4.4) 


Substituting this equation in the first two equations 
of equations (2.4.3), one gets. 


-(CA13-CA33CA93 CA73 ) - ( CAE3-CA33CA9T '^CASS) 
(CA43 -CA63CA93“^[:A73) ( CA53-CA63CA9n~'^ CAST ) 


0 r 


< 



1 

y 

1 



[ ^'^^2 


■■ CBIT 

CB2T ■ 


f 

- 




' ■fq>^ ' 

CB3T 

CB4T 


< 

CTJ^ 

y 

““ 

-< 

■Cq>2 


(2.4.5) 


(2.4.6) 


These equations are rearranged as 
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^q>. 


■CBZl ^cB-n 


CBZl 


-1 


CB33-CB4DCB23 ‘’cB13 CB43[:B2D 


•Cq > 


o r 








^ = CTFM3, . 


<q>2 

> 

1 



(2.4.7) 


Case Cb5 ; Convection on one edge and temperature prescribed 
on the other edge 

5 

4 
3 
2 

1 

Fig. 2. 4. 2 : Plane lamina with convection on one edge and 
temperature prescribed on other edge 




FEM equations for the strip 1 are the equations (2.4.1). 

To convert these equations into TFEM equations, the boundary 
conditions on bottom edge 1-11-6 (i.e. ^ - T^, - 

T. and B* = f? ~ T ) are applied and applying convection boundary 

T O O T 

condition on the top edge 5-13-10, one gets 
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Transfer matrix for this strip can be obtained from second, 
third, fourth, fifth, seventh, eighth, ninth, tenth, twelfth and 


E8 



thirteenth rows. Any one of the dropped out rows can later be 


used to find the flux along edge 16. 
mentioned ten equations 


' 22 ' 


' 32 ' 


' 42 ' 


' 52 ' 


' 72 ' 


' 82 ' 


' 92 ' 


' 10 , 2 ' 


' 12 , 2 " 

' 13 , 2 * 

10x11 


Thus writing the above 




e. 


9 


2 , 13 

'"21 ^1 " 


^ 2 , 11^11 ^ 


3,13 

<^31 < ^ 


S , 11 ^ 11 > 


4 , 13 

'Si ^1 " 

UbK * 

S , 11 ^ 11 > 


5 , 13 

'Si ^1 " 

a. 

S , ii < i ' 


7 , 13 

'Si ^1 ^ 


S , 11 ^ 11 ' 

< 

8 , 13 

'Si < " 


S , 11 < 1 ^ 


9 , 13 

'S i^i ^ 

* 

^9 , 1 1 *^ 1 1 ^ 


10,13 

' So , i <" 


e '' ) 

10,1 111 


12,13 

' S 2 , 1 <" 


9^^ ) 
12,11 11 


13,13 

' S 3 , 1 <" 


) 

13,1 r 11 _ 

L. 

r 

-^^31 

1 





(2.4.9) 
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Temperatures at the intermediate nodes of the strip are not 
required in the TFEM equations of the strip. These can be 
condensed out, C19D. 

Equations (2.4.9) are partitioned and rewritten as. 


CA11 

CA23 

CA30 

CX13 

4;<4 

4x4 

4x2 

4x1 

[:A40 

CA50 

CA63 

•CX2} 

4x4 

4x4 

4x2 

4x1 

CA70 

CA80 

CA93 

(1X3} 

2x4 

2x4 

2x2 

2x1 


iey 


Cd> 


-tq^ 


•Cd>, 




•CO} 


(2.4.10) 


Third equation of equations (2.4.10) gives, 
1 


Cd}^ = CA90 


j^-CA73 C^}.^ - CA80 •Cd}^ - CXS} j (2.4.11) 
Substituting this equation in the first two equations of 
equations (2.4.10), one gets 

-(CA10-CA30CA90“'^CA70) - ( CA20-CA30CA90 CASO) - ( ■CX1}-CA33CA93 ^ -CX 

(I:A40-CA61CA90''^CA73) ( CA5n-CA63CA93“ ^CASl ) ( •CX2}-CA8nCA93 ^ ■CX3}) 


o r 


(2.4. 12) 


(2.4.13) 




f r -k ^ 

■ce>2 

. = 1 


1 

V > 




■ CB10 

CB20 

■CB3} ■ 

-1 

f 

ret 1 'I 

'■®^1 

f 



f r 

CB40 

CB50 

■CB6} 







•Cq >2 



- 


< 





* 
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These equations are rearranged as. 


1 

r 

2 



ii 

^q>2 

V. J 

L 


-E:B2D ^CB13 


[:B2D 


CB23 iBSy 


(CB4D-CB53CB2D^ CB1D) I;B5DCB2D’^ ( CB6D-CB5DCB23 ^ CB 




iqy 


1 

1 y 


(2.4.14) 


Making state vectors (on left and right sides of the strip) 
of the same type, these equations get modified as 


•c©>„ I 

r 

2 



It 

■Cq >2 

V , J 

L 


-CB23 


LBZl 


(CB4D-CB53CB2D EIB 1 3 ) CBSDCBBD ^ 


LOJ 


LOJ 


CB23 '^CB3> 

( CB6D-CB53CB23“’^ CBSD ) 
1 


r 





0 r 


(2.4.15) 


1 


•cCT j- 

2 


1 

•C q > 2 

. = CTFMl,^ ^ 


1 

V. J 


1 

V. y 


a. 5 TFEM equations for stress analysis 

In this section TFEM equations are developed first for a 
domain with four noded finite elements and then with eight noded 
finite elements. 
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CAseCaUi Domain with four noded flnito elements 

YA 


x,u 

Fig. 2. 5. 1 : Rectangular lamina with pressure p applied on its 
upper edge 

As before, lamina is first divided into number of strips. 
Strip can be divided into any number of finite elements. For 
reducing typing work and easy understanding, each strip is 
supposed to have two elements only. 



FEM equations for a four noded element are C21D 


(ne) _ j.pj(ne) ^ ^p.,(ne) 


(2.S. 1 ) 


where , 

is the elemental stiffness matrix 
^pj<ne) elemental load matrix 

^Fj<ne) force vector along the edges of the element 

FEii equations for the first strip with two elements become 
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’ ^11 

K 1 2 • ■ ■ 

^ 1,12 


^1 

^21 

^ 22 * ■ • 

Kg , ^2 


'^l 

^31 

K32... 

^ 3,12 


^2 

^41 

K42. . . 

K4, ^2 


^^2 



m 

< 

U3 1 

''3 





^4 




! 











^5 



■ 


^6 

. ''■ 12 , 

1 

■ • ■ K .^ 2 , 12 _ 


L 


0 


■ ^^14 ^=^21 ■ 

0 


Fy,4 - Fy^.^ 

m 


FX 32 + Fxg.^ 

m 


*^^32 ^^21 

m 


^=^32 ^=^63 

< . 

> < 
+ 

■ ^^32 ^^63 

^'^14 *^'^45 

■ 


f"yi4 ^"^45 

li 


^^45 ^=^56 

11 


^^45 ^^56 

M 


^=^56 ■" ^^63 

m 

0 _ 


■"^56 '"^63 

L 


( 2 . 5 . 2 ) 


Since TFEM equations are to be developed from equation 
( 2 . 5 . 2 ), the edge forces ^^ 45 '^'^ 56 ' ^^56 changed to 

- Fx ,,., - FVc :., - Fx,cr ~^ y / i =* Noting the boundary conditions on 
34 I>4 03 o3 

the bottom edge 14 ~ '^°P 36 

0 ; f ^ y3 £, t»y Py3 ^^6^' equations ( 2 . 5 . 2 ) 
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^11 

Ki 2 . 


S 6 

^17' ■ 

• ’^ 1,12 

0 ■ 


■ ^1 ■ 



^21 

S 2 • 

••■S5 

S 6 

S7-- 

• *^ 2,12 

0 


'^i 


'^^21 

. 

• 





• 


■ 


+ Fx« . 

•» 






- 


m 


32 21 

Si 

S 2 • 


S 6 

S7-- 

■•S ,12 

0 


^3 


^^32 *^^21 











F x^« 

^61 

S 2 • 

. . 


S7-- 

■ ‘^ 6,12 


< 


>”< 

32 

*^^32 

si’ 

^72 ■ * 


S 6 

K77. . 

■ *^7,12 

0 


!■”■ 


-FK 54 

Sl 

^82' • 

• ■ ‘Ss 


S7‘ ■ 

■ ’^,12 

0 


% 


-^S4 

* 






• 


m 


-Fx^ .-Fx , 

• 






m 


m 


54 ^65 

11,1 

Sl, 2 - 

•Sl,5^ 

11,6 

^12,7" 

■ '^ 12,12 

m 

> 0 


m 

^6 


“‘^S4~'^^65 









V, 


■^S5 

12 , 1 

^ 12 , 2 ‘ 

•^12,5^ 

12,6 

^12,7- 

■^ 12, 12 



0 

1 


~^^65 


(2.5.3) 


Partitioning equation (2.5.3) as indicated, one gets. 


■ CAID i:A23 CXIJ ■ 


{u> . 


<F>, 1 



1 


1 


«< 


h = H 

► 

CA33 CA43 tX2> 


•Cu}^ 





a 


Z) 


1 


(2.5.4) 

As explained earlier, these equations can be rearranged as 





c = 

CF>2 


k. ^ 

1 



-CA23 ^CA-n 

(CA3D-CA4DCA23 

LOJ 


CA2D 

CA1D)-CA4DCA23 

LOJ 


-CA2n ^CX1} 

( CX23-[;A4nCA2D~'^ C X 13 ) 




K r -1 

tu>^ 




tF>2 

^ = CTFM3.^^ 


^ 1 


1 J 


(2.5.5) 


34 



Equation (2.5.5) gives the relation between the 
displacements and forces on the right side of the strip and the 
displacement and forces on the left hand side of the strip. 

CaseCbD : Domain with eight noded elements 



Fig. 2. S.a I Rectangular lamina with P o*' 

upper edge 

The finite element equations for the fii"st strip are 
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(2.5.6) 

Noting the boundary conditions on the bottom edge 16 (Fi<../ 

I O 

= Fy,^^ = 0)and on the top edge 5,10 (F;<^ = 0; Fy^ by Py^ 

Pyi 2 3ind Fy^Q^ » equations (2.5.6) become 
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‘‘•11 

*^12 

*" 1,26 

0 


“l 


'">‘31 

‘‘•21 

‘‘•22 

’' 2,26 

0 


'^1 



‘‘^31 

‘^32 

‘‘• 3,26 

0 


^2 


'">‘31 

‘‘^41 

‘‘^42 ••••■■ 

‘‘• 4,26 

0 




‘"ysi 

‘'^51 

•<52 

‘" 5,26 

0 


^3 


"^^53 "■ '^="31 

‘^61 

‘‘^62 

‘‘• 6,26 

0 


"-3 


^^53 " ^^731 

‘‘ T1 

‘^72 

‘" 7,26 

0 


^4 


'"“53 

‘‘•81 

‘‘■82 

‘" 8,26 

0 


'^4 


‘^^53 

‘‘■91 

‘‘^92 

‘‘• 9,26 

0 


""5 


FX53 

L : 

10,1 

‘‘^ 10,2 

‘" l0,26 

- Py5 


"^5 



‘'■ 11,1 

‘" 11,2 

‘‘H 1 , 26 

0 


^6 


-'">‘86 

U - 

12,1 

‘‘• 12,2 

‘‘• 12,26 

0 


"'6 


“*^^86 

k 

k - - ^ - 

k 

0 

< 

Q 

^ < 

— F K ,.4 , 

13,1 

13,2 

13,26 



7 


86 

‘‘• 14,1 

*' 14,E 

‘" l4,26 

0 


^'7 


-^786 

‘‘■ 15,1 

‘‘^ 15,2 

‘" 15,26 

0 


'•‘s 


-^786 - ‘^" 10,81 

‘‘■ 16,1 

‘‘■ 16,2 

‘" l6,26 

0 


'^8 


'^^86 ~ ’^ 7 . 10,8 

‘‘^ 17,1 

‘‘^ 17,2 

‘" l7,26 

0 


•^9 


“■‘^^ 10,8 

‘‘■• 18,1 

‘‘^ 18,2 

‘" is, 26 

0 


'^9 


”‘^ 7 . 10^8 

‘‘• 19.1 

‘" 19,2 

‘‘■ 19,26 

0 


^10 


”'^^ 10,8 

‘‘^ 20 , 1 

‘‘^ 20,2 

‘‘• 20,26 



"'lO 


“‘^ 7 ^ 0,8 

t- 


k . 

0 


U , , 


0 

21 , 1 

‘‘^ 21,2 

21,26 



1 1 



‘‘• 22, 1 

‘" 22,2 

‘" 22,26 

0 


'^11 


0 

‘‘• 23,1 

*" 23,2 

‘" 23,26 

0 


^12 


0 

‘‘^ 24 , 1 

‘" 24,2 

‘" 24,26 

0 


'^12 


0 

‘‘^ 25 , 1 

*" 25,2 

‘" 25,26 

0 


“13 


0 

‘‘• 26,1 

*" 26,2 

‘" 26,26 

-py-13 


, ''is 

1 


0 






L ^ J 


( 2 . 5 . 7 ) 
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are 


not 


Since displacement at the intermediate nodes 11,12,13 
required in the TFEM equations of the strip, these can be 
condensed out C19I]. Thus partitioning the equations (2.5.7), and 
rewriting them , 


CA 1 3 

10;<10 

CA23 

10x10 

CA33 

10x6 

CX13 

10x1 


■ -'"'l ■ 

10x1 


■ ■ 

10x1 

CA43 

10x10 

CA53 

10x10 

CA63 

10x6 

-CX2J 

10x1 

j 

•Cu3g J 
10x1 


-CF>2 , 
10x1 

CA73 

6x10 

CA83 

6x10 

CA93 

6x6 

•CX3> 

6x1 


6x 1 

1 J 


coy 

6x1 


(2.5.8) 

Third equation of equations (2.5.8) gives terms of 

■Cu>^,and £u}^. These can be substituted in first two equations 
of (2.5.8) , one gets 

(CA13-CA33CA93“"'CA71) ( I:AE3-CA33CA93'^CA83 ) ( CX 1 >-CA33CA93“ CXSJ ) 
-(CA43-CA63CA93“’^CA73)-(CA53-CA63CA93“'^CA83)-( fX2>-CA63CA92''’^-CX3> ) 




-F-l 

- -fUj* 


)<• 

2 

1 


fF^ 

^ * J 




(2.5.9) 


or 


CB13 CB23 
CB43 CB53 
These equations are similar to 
rearranged in the same way as. 


•CB3> "j r tuJ. T 
\ CuJ [ 

fB6> J t ^ J 



(2.5.10) 


equations (2.5.9) and are 
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' -CuJg ' 


IFJ^ 

K — 

2 


1 



■ -CBH3“''CB13 
( CB4D-CB5DCB23 
LOJ 


CB1D) 


-CB2n“''^B3> 

CB53CB2D"'^ (CB6D-CB5DCB23“'’-tB3 

LOJ 1 


Cii> 

■CF> 


(2.5. 11) 


0 r 




•Cu>.^ 

■CF>2 

1 

. = CTFMJ^ 

1 

CF> . - 

1 

1 


(2.5. 12) 


These are the transfer FEM equations for the strip 1 


2.6 TFEM for heat conduction problems 


y4 


Fig. 2. 6.1: Rectangular lamina divided into n-strips 
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In articles 2.3 and 2.4 transfer finite element 
were developed for certain heat conduction problems 
typical strip, oine get 


equat i ons 
For a 
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•CT>2 




' = CTFM3 , 

► 

1 J 

1 



( 2 . 6 . 1 ) 


Referring to Fig. 2.6.1, for a domain with n strips, on repeated 
multiplications one reaches the last strip as. 


■CT} 


LqJ- 


n+1 

n+1 


= CTFM3 CTFMD , 
n n-1 


CTFMH 


•tT> 
•t q } 


o r 


r 1 

n+ 1 


1 


' -CT} 

•CqJ 

n+1 

- = CTFM3„. , . 

final 

•tq3-i 

► = Lzl ^ 


1 


1 

k J 


1 

k J 


( 2 . 6 . 2 ) 


Here CzH is the final transfer matrix. Now the boundary- 
conditions on the left side of the first strip an d on the right 
side of the last strip are applied. Thus one can find the 
temperatures and heat fluxes. This is explained here for one set 
of boundary conditions, and extension to other sets of boundary 
conditions will be obvious. Here case (d) of article 2.3.3 is 
taken i.e. fluid on the upper edge and prescribed temperatures on 
the lower edge. Its TFEM equations will be similar to those of 
equation (2.3.26). Referring to Fig. 2.6.2 let the boundary 
conditions on the left hand edge the prescribed temperatures and 
on the right hand edge be insulation. 
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Tfluid 



Fig. 2.6.2 


Rectangular lamina with boundary conditions 
prescribed on the first and the last edge. 


Considering -Pour strips in the domain the final TFEli 
equations (2.6.2) can be written as 


' ^14 ■ 


*y 

11 

^12 

^13 

^14 

^15 


■ ^2 ' 

^15 


21 

^22 

^23 

^24 

^25 


^3 

■^^ 115,13 


*r 

31 

^32 

^33 

^34 

^35 

< 

■^^Isi j> 

15,14 


^41 

^42 

^43 

^44 

^45 


q 1 1 32 

■ 1 


1 

0 

1 

0 

0 

0 

1 


1 

> 


L 






■i Ml 


(2.6.3) 

Boundary conditions are Tg = Tg, Tg = Tg and ql|^^ ^g = 0, 

= 0- With these third and fourth equations give 
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r 



^33 

’34 

-1 

-v 

— z T — T T — 

31 2 32 3 35 


qi j 32 

K — 

*•43 

^44_ 


T* , T* - 

^41 2 ^42 3 ^45 


(H.6.4) 

After having found the flux values qi|3i qM 22 ' ^ tate 
vector on the left side of the first strip is completely known. 
Thus using the TFEM equations of each strip, starting with the 
first, one can find temperatures and fluxes everywhere. 


2.7 TFEM for plane stress problems 
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Fig. 2.7.1 : Rectangular lamina divided into n-strips and 

stibjected to load p 

In article Z.5 TFEM equations were developed for plane 
stress- problems, see equation (2.5.5). 

Referring to Fig. 2.7.1, for a domain with n-strips, on 
repeated multiplication one reaches the last strip as 





■Cu> 


. = CTFMT CTFMT^ . CTFMT ■ 

n n~ 1 1 


CF> ' 

1 

1 1 J 
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or 



Now the boundary conditions on the left side of the first 
strip and on the right side of the last strip are applied. Thus 
one can find the forces and the displacements. This is explained 
here for one set of boundary conditions, and extension to other 
sets of boundary conditions will be obvious. 

Here case (a) of article 2.5 is taken. Let the first and 


the last edges be fixed in both directions. Fig. 2.7.2. 

P 


y 

X 

Fig. 2.7.2 : Rectangtilar plane wibh the first and the last edge 

fixed. 

Taking three strips in the domain, the final transfer matrix 
eg uat ion (2.7.2) can be written as 
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^10 


"11 "12 ^1,13 



''10 


"21 ^22 -2,13 


''i 

^11 


"31 -32 -3,13 


^2 

''11 


"41 ^42 -4,13 


"2 

Ui2 


"51 ^52 5,13 



''l2 

!> 

"61 -62 -6,13 


''3 i 

*^^11 , 10 


-71 -72 -7,13 


Fx2 .^ 

*^^11 , 10 


-81 -82 -8,13 


^^21 

’^'"'12, 10 


-91 -92 -9,13 


f^^<31 

*^^12,10 


-10,1 -10,2 ^10,13 


*^^31 

^'"'12, 11 


^11,1 -11,2 -11,13 


FX32 

^^12,11 


-12,1 -12,2 -12,13 


^^32 

1 


L 0 0 1 


1 


(2.7.3) 


Boundary conditions are u^ = v 


u 


2 C. O O 

With these first six 


and 


‘^■10 = ^^10 


^11 


^■12 = ''12 = ^■ 


equations give 


1 , 13 


-17 

-18 

"19 

2, 13 


^27 

"28 

"29 

3,13 


Z37 

^38 

^39 

4,13 


^47 

^48 

2 49 




“7 


5,13 


•^57 

"58 

"59 

6,13 


^67 

^68 

^69 


^1 , 10 

^1,11 

^1,12 


* Fx^.^ 

-2,10 

^2,11 

^2,12 


*^^21 



z 


Fx„ . 

^3,10 

"3,11 

3,12 

< 

^“31 

0 

N 

^4,11 

^4,12 

*^^31 

^5, 10 

^5,11 

^5,12 


FX32 

^6, 10 

^6,11 

^6, 12 


1 

OJ 

CO 

u. 

1 


( 2 . 7 . 4 ) 
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i.e. force vector on left side the first strip becomes known. 
Thus state vector on the left side of the first strip is 


completely known. Thus using the TFEM equations of each 
starting with the first, one can find displacements and 
everywhere . 


strip, 

forces 
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CHAPTER -3 

THERMAL AND STRESS ANALYSIS OF AXI-SYMMETRIC PROBLEMS 


Axi-symmet r i c bodies subjected to axi-symmetri c loading 
often have complex geometrical features. For example a vertical 
pressure vessel with a spherical or ,hemi -spherical or 
,tori-spherical or an elliptical head subjected to 
ax i-symmet r i c thermal and mechanical loadings. At times such 
vessels have openings (inlet and outlet) in the centre of the 
head and thicknesses of the head and shell are different. 

In this chapter TFEM equations will be developed for steady 
and transient heat conduction and stress analysis problems for 
axisymmetric bodies. Such problems are mathematically two 
dimensional problems. 

3.1 Governing differential equation 

The governing differential equation for heat conduction for 
axisymmetric problems are LZ31 


S 

Sr 


(k 




W 


ST 

-(k ^) = 

: z Sz 


P c 


ST 

St 


~ ( 3 . -1 . 1 ) 


with boundary conditions 


T=T (r,z) at t = 0 

o 


and 


k r 
r 


ST 

Sr^ 


Sz 


n 


+ q + h 


(T - ^ 


0 


(3.1.2) 


3.2 FEM equations, heat conduction 

Denoting the temperature difference T - by S, the 

governing differential equation (3.1.1) can be rewritten as. 


46 



r r ar Sz^^z Sz^ ~ ^ ^ 5t 

and boundary conditions become 


(3.2.1 ) 


^ ~ t = 0 


and 


. se ^ , ae 

57 5z "z 


n_^ + q + h& = 0 


(3.2.2) 


(e) 


Assuming d over the element 

= LN(r,z)Jtd(t)>^^®^ (3.2.3) 

and substituting in equation (3.2.1), one gets the residue as 


( e ) 


i £_ (I. r ^ 

r ar r Sr 


e) 


^(k 

oz z 


SB 

Sz 


( e ) 


. SB 


(e) 


(3.2.4) 


Minimizing this residue by Galerkin method, one gets 


r r ht r ’1 ^ /I ^ ^ \ u. 

II N . - -3— ( k . r . -3— ) + -3— 

JJ ilr^r r Sr Sz 


^ (e) 

^^•z 57 I'' dr dz=0 


Doing integration by parts, one gets. 


j; [<" 


. r . CN , r } -r- + k . r . CM , z 3-^ 

r Sr z Sz 


) — p cr ■CN> 


(e) 


se 

St 


(e). 


^ - CN3 ( q + h d "" ' ) dB 


dr dz 


(3.2.5) 


Substituting equation (3.2.3), these become 
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.r.-CN,r>LN,rJ + k r.-CN, 2 > lN,zJ ) -C^> 


IS [‘^ 

+ p cr -CNJ LNJ ■Cd>^"®^jclr dz = -CNJ <q + h lNJ d ^ ) dB 


(3.2.6) 


Equation (3.2.6) in short form is. 


CCl 


(e ) 


■W ) 


( ne ) 


Ck^n 


(e) 


•Cd} 


( ne ) 


Ck. □ 
h 


(e) 


■ze: 


(ne) 


= CqJ 


( ne ) 


(3.2.7) 


where 
( e ) 


CCD 


CktD 


(e) 


Ck. D 
h 


( e ) 


- , (ne ) 
■Cq> 


= j jp cr •CN> LNJ dr dz= Thermal capacitance matrix 
= J J (k^.r.-£N,r} LN,rJ + k^ . r . CN , LN , ^ J ) dr dz 
= Thermal stiffness matrix 

JL hCN> LNJ dB = Elemental convection mofcAtx where 

appl i cable 

= ^- q.CNl dB = Elemental flux vector where applicable 


For steady state heat conduction problems, equation (3.2.7) 
becomes , 

Ck,D <3.2.8) 

tr 

3.3 FEM eqtiationst, ssbress analysis 

For stress analysis, variational approach is more convenient 
and is used here. The displacements in an eight noded element 
case are approximated over the eight noded element as. 
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N.^ 0 
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•f u> 


( ne ) 


CND 


•lUw 


( ne ) 


The strain matrix is C22I1 


(3.3. 
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► 1 



^v/Sz j 
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^ STr J 


Wr~ 




0 

C. 
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0 

CL 

r 

0 

0 

#<^•1 

Sz 

0 
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sT~ 

dr 

az 

dr 


8 ^ 1 


r 1 



u , 

dr 


1 

No 


V ^ 

8 y-V 

2 0 


1 

r 


^2 

dNo 



0 , ® 

< 

2 

dz 


m 

m 



m 

m 

^8 

dz dr 


'^8 


(3.3.2) 


and the stress matrix is written as 


1) 
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CDD L^l 


(3.3.3) 


i&y = 

where CD3 is the elasticity matrix C213. 

Using the above equations, the potential energy for 
ax i-symme t r i c elements becomes 


,,(e) 1 , ,(ne) 

W = = LuJ 


J J CBl''" CD3 CB3 2 fi r dr dz 
( e ) 


A 


A<e) 


2 n r drdz - JtNl 


T 


V. J 


2 n r dB 


(3.3.4 


where F^, are body forces per unit volume and q^, are 
surface forces per area. Minimizing this potential energy one 
gets 


J J CBl'’’ CDD CBl 2n r dr dz ■Cu>^'^®^=J J Cn 5 
(e) 


A 


A 




V J 


Z n r dr dz + J CNl 


{:: )'■ 


rdB (3.3.5) 


and this in short form is 


50 



Ck3^®^ (3.3.6) 

B s 

where 

( 0 ) 

Ck3 is the elemental stiffness matrix 
•CPg>^'^®^ is the elemental body force vector 

■CP is the elemental surface (external or internal 

s 

force vector) 


3. it TFEM equations for 'the steady state heat conduction problesas 

Equations (3.2.8) for this heat conduction are similar to 


equat i on 

(2.3.6) . 

Thus TFEM will 

be 

simi lar 

to equations of 

article 

2.5 for 

various types 

of 

boundary 

conditions . 


3.5 TFEM equations for transient heat conduction problems 



r 


Fige3. 5 : Axi-syiranetrie strip with convection on the inside edge 

and temperature prescribed on the outer edge 

113493. 
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Figure 3.5 shows a typical strip as used in the 
ax i-symmet r i c problems, composed of two ring elements with fluid 
on the inside edge and prescribed temperature on the outer edge. 

To deal with this transient heat conduction problem governed 
by equation (3.2.7), the use is made of the most successful 
scheme of 'Implicit Scheme’, thus -tdl is represented as 


. ( ne ) 

i&y 


) 


(ne ) 

p+1 


■Cd> 


(ne > 
P 


At 


(3.5. 1 ) 


where p is the previous time 
p+1 is the current time 
At is the time step 
Equations (3.2.7) became 


( , cc n 

AX 


+ + 


(k,3) 

n 




p+1 


— .r .-I \ 4- 

' ""'^p+1 


1 

At 


cci] wy 


(3.5.2) 


Using the Implicit scheme illustrated by equation (3.5.2), 
and noting convection on the edge 5-13-10, the FEM equations are 

given by 
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(3.5.3) 


The state vector on the RHS in equation (3.5.3) is 

known from initial conditions. Denoting the matrix 

(74-i:C3+Ck,n+Ck^3) by Cr3, and the vector CCJ i&J by -CxJ, 

At th-^ At p ' 

the finite element equations (3.5.3) after applying boundary 
conditions of the outer edge and rearranging become 
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13x1 


TFEM equations for this strip will be developed from second, 
third, fourth, fifth, seventh, eighth, ninth, tenth, twelfth and 
thirteenth equations. Any one of the dropped out equations can 
later be used to find the flux along edge 1-11-6. 
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Rewriting these equations. 


• 10 , 2 ' 


''2,-l3 1-1'"''26® 6"='2^ 

■^3, -13 ^ ''S 

'^4, •13^''4-1^1'^'^4, -M^11’^''46^6~'^4^ 

5 , -1 3 [ 5 -1^ -1 5 , -M^ 11 ■*■ '^ 5 6^ 6 “ 'S ^ 

''y, 13^ ‘^7-1^ I'^'^y, 1-1^1 l‘^''76^I 
‘'a, 13^ '^s-i^i'^'^a, i-i^ii'^'"86^6"'^a^ 

'^9, 13^ ‘'9-1^ l"^*^?, 1-1^1l'^''96^6~'^9^ 
*^10, 13^ “"lO, l^l'^'^IO, 1-1^1l'^''l0,6'^6 
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, V.'* v.-**- 10' 
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*"13, 13^ ''13, 1^1l'^'"l3, 1-1^1l‘^''l3,6^6 

-=<13^ 


10x11 


11x1 



10x1 


(3.5.5) 
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These equations are similar to eqns 
developed in a similar manner. 


(2.5.9) 


and TFEM can be 



(3.5.6) 


Thus i&y and Cql can be determined at all nodes at current 
time. Thus procedure is repeated for the next time step and so 
on . 


3. 6 TFEM equations for stress analysis 



r 

Fig. 3.6 ; Axisymmetrlc strip with loading p on the inside edge 


Figure 3.6 shows an axisymmetric strip subjected to a load p 
on its inside edge. The formulation for this case is quite 
similar to the case (b) of article 2.6. Rewriting the FEM 
equations for this strip after applying the boundary conditions 
along the edge 1-11-6 ~ *^^16 ~ edge 5-13-10 
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Since displacements at the intermediate 


nodes (11,12,13) 


are not required in the TKEM equations of the strip 
condensed out. 

Partitioning equations (3.6.1) as indicated 


them. 


CA13 

CA23 

CA30 

■CX1> 

10x10 

o 

6 

10x6 

10x1 


CA53 

CA6n 

•CX2> 

o 

o 

o 

6 

10x6 

10x1 

CA7n 

CA83 

CA93 

■CX33 

6 X 1 0 

6x10 

6x6 

6x1 


these can be 


and rewriting 




•CF>., ■ 
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\ ^ ( 

> < 

2*' 

10x1 

— 

10x1 



■COJ 

6x1 


6x1. 

1 




(3.6.2) 


After condensing the first two equations of equation 

(3.6.2) can be written as. 


■ CB13 

CB2D 

CB33 ■ 


1 1 


<F>^ ■ 

10x10 

10x10 

10x1 
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CB5n 

CB6n 


•Cu>2 


•CF>2 

10x10 

10x10 

10x1 


10x1 

1 


10x1 


(3.6.3) 
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One finally gets the transfer matrix as. 
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3. 7 


a general case 


TFEM equation for heat conduction for 

Z ^ 


© 


Tfluid 

@ 



Fig. 3. 7 : A general axi symmetric case 


Fig. 3.7 shows a general ax isymmet r i c case. Development of 
TFEM equations in such problems is different from the previous 
cases, owing to the presence of more elements in some strips than 
other strips. Here strip E is shown to have two additional 
elements in the z-direction. The additional nodes in such strips 
are condensed out. 

TFEM equations for strip 1 after applying boundary- 
conditions along edge 14 and edge 36 are 


r 1 


r 1 

4 




= CTFMl.^ 

^2 1 

] 1 

> < 

qil .|3 

‘^^l46 


13 

1 


1 
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Finite element equations of strip 


after applying the 


boundary conditions of edges 47, 6-27, ZJ-ZB and 28-9, are 


^44 ^45' '^4, ■19'‘^6^46^21^4,21'^23^4,23'^’ ‘ ■'*^9^4,9' 

^54 *^5, 19^^6^56‘^21^5,21'^23^5,23'*‘’ ■ •‘^9^5,9’ 

^74 ^75- -■^7, ■19^^6V6'^21 V,2-1'^23V,23‘'- • •’^9^7,9^ 
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7x1 


Since strip 2 contains 


(3.7.2) 

additional intermediate nodes. 


so 
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they will be condensed out to make TFEM equations compatible for 

marching. Thus condensing nodes 13,16 and 19, the TFEM equations 
for strip 2 are : 


■ ■ 

<, '^^79 

[ = CTFMD 

■ ■ 

1 

■^M79 

1 


1 


Similarly TFEM equations for strip 3 are given as. 


■ ^io' 


^7 



& 

1 1 


a 

^ 10-12 

> < 
= CTFMl^ 

•31 1 7, 

10-1Z 


■’477 

1 


1 


Now transfer matrices of strips i.e. equations (3.7.1), 
(3.7.3) and (3.7.4) can be multiplied and results obtained as 
usual . 


3.8 TFEM for axisymmebric problems 

TFEM equations for a typical strip have been derived in 
articles 3.4, 3.5, 3.6 and 3.7 for axi-symmetr i c heat conduction 
and stress analysis problems. Now one can multiply the transfer 
matrices for various strips and get all the results after 
applying the boundary conditions on the first and last edge, as 
explained in Chapter Z. 
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CHAPTER - 4- 


RESULTS AND DISCUSSIONS 

Computer programs in FORTRAN language have been made 
■for plane lamina and axi — symmetri c thermal and stress analysis 
problems incorporating different boudary conditions. A 

rectangular domain is considered in case of plane lamina piroblems 
and a pressure vessel with hemi— s phe r i cal heads is considered 
is case of ax i-symme t r i c problems. 

4,-1 Thermal analysis of rectangular lamina with four noded 
rectangular elements 

To validate the TFEM, first results were obtained for 
rectangular laminas with four noded rectangular elements. Each 
lamina was 5x'10 m. Coefficient of conductivity was taken as 54 
W/m. C. Each lamina was divided into ten strips and each strip 
in five elements, as shown in Fi gs . 4. -1 . '1 , 4. -1.2, 4. ‘1.3 and 4.1.4. 
As said in Chapter 2, following boundary conditions are 
considered on top and bottom edges. 

(a) Top and bottom edges insulated. 

(b) Top edge with prescribed temperature (200°) and bottom 
edge insulated. 

(c) Top and bottom edges with prescribed temperatures i.e. 
300°C at top and 100°C at bottom. 

(d) Top edge exposed to fluid (convection) at 100° and 
bottom edge with prescribed temperature (•10°C). 

Following boundary conditions are considered on the first and 
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Isst, edges for the above four cases. 


(a) First edge with prescribed temperature SOO^C and 
last edge with prescribed temperature of •10°C. 

(b) First edge insulated and last edge with prescribed 
temperature of lOO^C. 

(c) First and last edge with prescribed temperature of 
•100°C. 


(d) First and last edge insulated. 

Figure 4.1.1 shows results for case (a). Temperatures along 
the vertical edges of strip are identical, as this problem is 
essentially one dimensional heat flow in x-direction. These 
results obtained by TFEM match exactly with the analytical 
results C203. 


Results for case (b) are shown in Fig. 4.1.2. This is a two 
dimensional problem. Temperature decrease from left to right and 
also from top to bottom, as expected. Results for case (c) are 
shown in Fig. 4.1.3. This is also a two dimensional problem. 
Temperatures are symmetric about the middle vertical edge, and 
are maximum along this edge and decrease on both sides (left and 
right); as temperatures are same on the first and last edge. 
These temperatures decrease from top to bottom as expected. 


Temperatures for case (d) are shown in Fig. 4.1.4. 


Heat 


2 o 

transfer coefficient for this case is taken as 50 W/m . C. 
Temperatures along the horizontal edges of the strip are 
identical, as the problem is essentially one dimensional heat 
flow in y-direction. These results fay TFEM match exactly with 
the analytical results C203. 
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Fig. 4.1.2 Top edge with prescribed temperature and bottom 
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Fig. 4. 
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prescribed temperature. 


4.2 Thermal steady and transient analysis of pressure 


vessel with eight noded elements 

Vertical pressure vessel with hemispherical heads shown in 
Figure 4.2.1 is studied. Inside diameter is 2000 mm, thickness 
is 100 mm and outside length is 4200 mm. 

First following two cases of steady state heat conduction 
are studied. 

(i) Temperature prescribed at inner surface as 100°C and at 
outer surface as 10°C. 

<ii) Inside surface exposed to fluid at 250 C and outside 
surface with a prescribed temperature of 30*^0. 

Coefficient of conductivity was taken as 54 W/m. C and 
coefficient of convective heat transfer as 5000 W/m . C. It was 
noted that because of symmetry there will be no heat flow across 
the lines r and r^, i-e. q = 0. These are the boundary 

conditions for the first and last edge. 

Vessel was divided into 10 strips with two eight noded 
elements in each strip, as shown in Figure 4.2.1. 

Temperatures obtained by TFEM for case (i) are shown in Fig. 
4.2.1. These are symmetric with respect to the central 
longitudinal axis of the pressure vessel and also w.r.t line r^ as 
expected. Temperatures at node points along the line r^ in the 
head by TFEM are 75.832, 52.845, 30.911 and by analytical method 
applied to sphere alone are 75.853, 52.857, 30.930. Temperature 
at node points along line ,3 in tho shell by TFEM are 76.650, 
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53.905, 31.679 and by analytical method applied to cylinder alone 
are 76.683, 53.938, 31.708. Also temperatures decrease along the 
radial line, as expected. 

Temperatures obtained by TFEfI for case <ii) are shown in 
Fig. 4.2.E. They follow the trend of case (i). Temperatures at 
inside surface of the pressure vessel are less than the fluid 
temperature as expected. Values at the nodes along the line r .j 
by TFEM are 2S7.279, 174.335, 123.997, 75.874 and by analytical 

method for sphere alone are 226.641, 173.834, 123.640, 75. 7-^2. 

Similar results for line r^ are 227.605, 176.340, 126.401, 

77.600, 227.608, 176.412, 126.450, 77.664. Thus it is concluded 

that TFEM results are quite accurate. 

Next case (ii) was studied again for a pressure vessel with 
skirt, as shown in Fig. 4.2.3. Here vessel is again divided into 
10 strips with the eight noded elements in each strip except in 
the strip 8 (having skirt) which has five eight noded elements. 
Results are shown in Fig. 4.2.3. These are quite close to that 
of the vessel without skirt. Fig. 4.2.2, except in the strip 8 
where results change quite a bit due to the effect of the skirt. 
Moreover the results are not symmetric with respect to the line 
r^, as vessel is no more symmetric with respect to this line. 

Next transient heat . conduct i on was studied in the pressure 
vessel of Fig. 4.2.1. Vessel has fluid at 250°C inside and its 
outside surfaoe is maintained at 30°C. Value of P is taken as 
7833 kg/m^ and specific heat C as 468 d/kg."c. Time steps 
considered were 5, 10. 20 and 30 seconds. Time steps of 5 and 10 


70 



seconds were leading to numerical instability after some time. 
Time steps 20 and 30 seconds were found satisfactory. Results 
are shown in Figure 4.2.4 at 20,60, 100, 300, 400 and 500 seconds 
using the time step of 20 seconds. They show the expected trend. 
At 500 seconds the temperatures are nearly the steady state 
values. 
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Fig. 4.2.1 Pressure vessel with prescribed tempera- 
tures at inside and out side surfaces. 
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Fig. 4 . 2.2 
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Fig. 4.2.3 Pressure vessel with skirt. 
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■4.3 Stress analysis pressure vessel with eight noded 


elements 

Vertical pressure vessel of Fig. 4.2.1 is studied again for 
stress analysis. Here only half of the vessel is considered 
because of the double symmetry. The same should have been done 
for the theriT.al analysis, but it got overlooked there. Inside 

pressure is 10 MPa. vessel was divided into 5 strips with two 

eight noded elements in each strip. 

Because of symmetry there will be no r-di s placements along 

line r^ and no z-dis placements along the line r 3 . This is shown 

in Fig. 4.2.5. Natural boundary conditions along line r^ are 

that forces in z-direction are zero, and along line r 3 forces in 

the r-direction are zero. These geometrical and natural boundary 

, ,t the end of the multiplications of 
conditions were used ai- ■' 

transfer matrices of all the strips. 

Displacements and stresses are obtained by TFEH. Radial 

displacements are shown in Fig- ''.E.i. They show the e«pected 

trend. Radial displacements by TFEM alon .3 the line r 3 are .4754. 

.4685, .4623, .4562, .4504 mm and by analytical method for a 

, ■ c, r ,, ro„-, 5411, .5362, .5316, .5275, .5238 mm. 

straight shell .C24D are 


These results do not ma 


tch very well as shell length is quite 


short and head geometry 


affects the displacements 


Vertical 


displacements are shown 


Fig. 4.2.7. These displacements also 


follow the expected trend. 

. the Gauss points. Radial stresses 

Stresses are obtainea 

t .idinal aauss lines in Fig. 4.2.8. 


, , innoitudinal gauss lin( 

are shown along the lonj^* 

t t rit-TTeenth ti'^nsverse ge-uss 

These values along fiTvet? 


line are ~‘9»814y 
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-7.195, -4.460, -4.597, -2.3E9, —.011 MPa and by analytical 
method for a straight shell CS43 are -9.354, -7.223, -5.210, 
-4.085, -2.240, -.492. These results are comparable. Radial 
stresses along the first transerse gauss line are of the order of 
50 MPa and compare well with the spherical shell theory value of 
50 MPa. Hoop stresses along the si:< longitudinal gauss lines 
are shown in Fig. 4.2.9. They show the expected trend. These 
values along the fifteenth transverse gauss line are 105.838, 
103.790, 101.987, 100.274, 98.505, 96.892 MPa and agree well with 


analy t i ca 1 method for a straight shell CE43 values of 104 . 592 , 
102.461, 100.448, 99.323, 97.478, 95.731 MPa. Hoop stresses 
along the first transverse gauss line are of the 
order of 50 MPa and by analytical method for these spherical 
shell are also of the order of 50 MPa. o'^-stresses along six 
longitudinal gauss lines are shown in Fig. 4.8.10. These also 
show the expected trend. These stresses along fifteenth 
transverse gauss line are of the order of 50 MPa and by 
approximate thin shell theory are also of the order of 50 MPa. 
Along first transverse gauss line, changes from -9.512 to 
-.069 MPa. This is also consistent as pressure inside is -10 MPa 


and 0 at outside. 
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Fig. 4.2.5 Pressure vessel subjected to an 
internal pressure. 
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Fig 4.2.10 Longitudinal stress vs Gauss points 


CHAPTER- 5 


CONCLUSIONS 

Following conclusions have been arrived at in the present 

work; 

i) The TFEM equations can be obtained for elements having any 
number of nodes and any degrees of freedom. 

ii) The transfer finite element method can be applied 
successfully to transient problems also. 

iii) The transfer finite element method can be applied to 
pressure vessels having any geometry as long as the vessel 
is ax i-symmet ri c . 

iv) Results obtained by TFEM are accurate enough to match with 


the exact results. 
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